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We compare phenomenological values of the frozen QCD running coupling constant (a s ) with two 
classes of solutions obtained through nonperturbative Schwinger-Dyson equations. We use these 
same solutions with frozen coupling constants as well as their respective nonperturbative gluon prop- 
agators to compute the QCD prediction for the asymptotic pion form factor. Agreement between 
theory and experiment on a s (0) and F- K {Q 2 ) is found only for one of the solutions Schwinger-Dyson 
equations. 



I. INTRODUCTION 

The possibility that the QCD coupling constant (a s ) 
has an infrared (IR) finite behavior has been extensively 
studied in recent years. There are theoretical arguments 
in favor of the coupling constant freezing at low momenta, 
one of them, a la Banks and Zaks [Ej, claims that QCD 
may have a non-trivial IR fixed point even for a small 
number of quarks (see, for instance, Ref. We can also 
use arguments of analyticity to show that the analytical 
coupling freezes at the value of 47r//3 || , where (3q is the 
one-loop coefficient of the QCD (3 function. 

Studies of the nonperturbative QCD vacuum also in- 
dicate the existence of a finite coupling constant in the 
IR |,|. 

The phenomenological evidences for the strong cou- 
pling constant freezing in the IR are much more numer- 
ous. Models where a static potential is used to compute 
the hadronic spectra make use of a frozen coupling con- 
stant at long distances 

Heavy quarkonia decays and total hadron-hadron cross 
sections arc influenced by the freezing of the coupling 
constant J§|. 

A quite detailed analysis of the ratio R e + e - (= 



c*ot(e + e — > hadrons) / <j(e + e 



)) performed by 



Mattingly and Stevenson [g| also shows a signal for the 
freezing of the QCD coupling. Following an almost sim- 
ilar study, for several hadronic observables, Dokshitzcr 
and Webber obtain the same result JlO| ]. 

Another method to investigate the infrared behavior of 
gluon and ghost propagators, and of the running coupling 
constant at low energies is through the solution of the 
Schwinger-Dyson equations (SDE) Early studies of 
the SDE for the gluon propagator in the Landau gauge 
concluded that the gluon propagator is highly singular 
in the infrared However, these results are discarded 
by simulations of QCD on the lattice at 95% confidence 
level [13], where it is shown that the gluon propagator 



probably is infrared finite. The lattice result is in agree- 
ment with two classes of SDE solutions. One proposed by 
Cornwall many years ago where the gluon acquires a dy- 
namical mass |14fl , and another that has been extensively 
discussed by Alkofer and von Smekal where the gluon 
pro pagator goes to zero when the momentum q 2 — > 
|15| . The solutions differ due to the different approxima- 
tions performed to solve the SDE, but in both cases there 
is a freezing of the coupling constant in the IR. 

Although the figures of the most recent lattice calcu- 
lation |13|] seems to indicate that the Cornwall's gluon 
propagator is the one that could better explain the re- 
sults, it is correct to say that the data is still not precise 
enough in the IR region to decide among the two possi- 
ble behaviors for the gluon propagators discussed in the 
previous paragraph. The purpose of our work is exactly 
to confront the IR values of the theoretical coupling con- 
stant, obtained with the solutions of the SDE, with the 
phenomenological data about the value of a s (0) in order 
to discriminate which one is the most suitable solution. 
Finally, these theoretical and phenomenological calcula- 
tions are outside the scope of standard perturbation the- 
ory, and a consistency check between them is the min- 
imum that we may require to know if these approaches 
make sense at all. In the next section we present the ex- 
pressions of the nonperturbative running coupling con- 
stant obtained with the SDE study, and compare them 
with some of the phenomenological values obtained for 
a s (0). In Section III we compute the pion form factor 
(Fit(Q 2 )) as a function of these coupling constants. It is 
known that F 7r (Q 2 ) is quite dependent on the behavior 
of a s at small momentum [ p^| . 

Therefore, this calculation provides a good test for the 
nonperturbative expression of the QCD coupling con- 
stant. Considering that solutions of SDE show a nonper- 
turbative behavior for the infrared coupling as well as for 
the gluon propagators, in Section IV we modify the ex- 
pression for the asymptotic pion form factor to take into 
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account these nonperturbative gluon propagators. In the 
last section we present our discussion and conclusions. 



II. THE PHENOMENOLOGICAL VALUE OF a s (0) 

At high energies it is believed that the property of 
asymptotic freedom allow us to perform reliable QCD cal- 
culations. However, the same is not true at low energies, 
where we have to make use of a series of phenomenolog- 
ical models when computing strong interaction parame- 
ters. 

This is exactly what happens if we want to determine 
the IR behavior of the running coupling constant. We 
are going to present some of the determinations of a s (0), 
and the most impressive fact is that the values obtained 
in several different analysis are not far apart by one order 
of magnitude, but they differ at most by a factor of two, 
providing a solid indication of the robustness of these 
approaches. 

One of the most detailed calculation of a s (0) is due 
to Mattingly and Stevenson || , which uses perturbation 
theory and renormalization group invariance to compute 
R e + e - up to third order in a s . They predict the value 



a s /ir = 0.26, 



(1) 



(a s (0) = 0.82) for the frozen IR coupling. On the other 
hand the long work of Ref. |n| gives 



a s (0) » 0.63. 



(2) 



The analysis of hadronic spectroscopy with potential 
models by Godfrey and Isgur H led to the following be- 
havior of the coupling constant 



a g i = 0.25 exp(— q 



0.15exp(-0.1g 2 
0.20exp(-0.001g 2 ), 



(3) 



where q is in GeV (all the momenta, otherwise specified, 
will be in Euclidean space), and a good fit of the spectra 
does not depend strongly on the ultraviolet behavior of 
the coupling constant. From the above equation we ob- 
tain a s (0) = 0.60. Which is also consistent with more 
recent studies of QCD potentials |l7| ]. Analysis of e + e~ 
annihilation, as well as bottomonium and charmonium 
fine structure in the framework of the background pertur- 
bation theory may lead to a frozen value of the coupling 
constant as low as a s (0) w 0.4 Jl^]. This method also 
explains the frozen value of a s resulting from the lattice 
simulation of the short range static potential |20| , and it 
gives 



a B {0) 



4tt 



A)ln- 



(4) 



where mg is a background mass. This one and Ay (with 
itib > Ay) are determined phenomenologically p0|. 



There are many other results that we could present 
here, but we can assume that the phenomenological val- 
ues of a s (0) scattered in the literature are in the range 



a s (0) « 0.7 ±0.3. 



(5) 



Although this choice is ad hoc, as far as we know it con- 
templates most of the phenomenological determinations 
ofa s (0). 

We now turn to the coupling constants obtained 
through the SDE solutions. The first nonperturbative 
running coupling constant that we shall discuss was ob- 
tained by Cornwall fbl , using the pinch technique to 
derive a gauge invariant SDE. This nonperturbative cou- 
pling is equal to 



a s c(q 2 ) 



-Itt 



/3 ln[(g 2 +4M fl V)/A 2 ]' 



(0) 



where M g (q 2 ) is a dynamical gluon mass given by, 



M 2 (q 2 ) = m 2 



In 



4m a ' 
A2 



-12/11 



(7) 



A(= Aqcd) is the QCD scale parameter, [3q = 11 — |n/, 
where n/ is the number of flavors. In the above expres- 
sion we are neglecting the effect of dynamical or bare 
fermions masses pfl . We can determine a s (0) in Eq.(j(]) 
as a function of the gluon mass m g and A, and these ones 
can be obtained in the calculation of several hadronic 
parameters that may vary with m g (but, in general, not 



strongly with the ratio m g /A). A typical value is 114 



500 ± 200 MeV 



(8) 



for A = 300 MeV. It is interesting to observe the similar- 
ity between Eq.(||) and Eq.(Q). Although, it is not clear 
to us the reason for this similarity. 

The other possibility for the IR finite running coupling 
was studied by Alkofer et al. ]l5| ], that solved a coupled 
set of SDE for the propagators of gluons and ghosts. In 
this approach the solution to the running coupling leads 
to an infrared fixed point, which, in terms of the invariant 
functions Z(k 2 ) and G(k 2 ) related to the renormalization 
of gluon and ghost propagators respectively, is given by 



Z( P 2 )G 2 (^)\^o 

-l 



4tt/3o 

16tt (1 1 
3N~ C Ik ~ 2 



9.5. 



(9) 



with k = 0.92. 

The above result gives a s near the origin. It has been 
obtained for rif = 0. As we are going to compare different 
SDE solutions we will limit ourselves to the flavorless, or 
pure gauge, QCD. The effect of n/ ^ will be discussed 
in the last section. 
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Since we shall use the running coupling in the full range 
of momenta we provide a fit for the numerical data of Ref. 
]l5| , given by the following expression: 



a s A 




q 2 < 0.31 GeV 

0.31 < q 2 < 1.3 GeV 2 

q 2 > 1.3 GeV 2 



(10) 



with 



a sA i = 0.2161 + 9.2621 exp 

a sA 2 = 1.4741 + 8.6072 exp 
1.4978 



(q 2 - 0.0297)^ 
(0.6846) 2 
q 2 - 0.1626 \ 
0.3197 J 



ln(1.8488g 2 ) ' 
where the x 2 ~ 2.5 x 10~ 4 for the three regions. 



(11) 




k [Gev*] 

FIG. 1. Comparison between the running couplings ob- 
tained from different approximations in the SDE study. The 
curves with line + triangle and line + circle delimit the phe- 
nomenological range acceptable for the gluon mass (300 MeV 
and 700 MeV, respectively) with A = 300 MeV for the Corn- 
wall's running coupling. The square points are the numerical 
data computed by Alkofer et al. and the solid line is our 
fit Eq.(^o[). The box in y axis shows the phenomenological 
range indicated in Eq.(^). We can see that only the running 
coupling computed by Cornwall is compatible with the phe- 
nomenological estimatives of a s (0). 

In Fig.(Q) we indicate the expected phenomenological 
range of values for a s (0) and plot the curves for a s c and 
a s A- 

It is evident that only the Cornwall's solution is com- 
patible with the phenomenological data. In the last sec- 
tion we shall comment on possible modifications of this 
result. 



III. THE NONPERTURBATIVE COUPLING AND 
THE PION FORM FACTOR 

It is known that the pion form factor, F^IQ 2 ), is quite 
dependent on the behavior of a s at small momentum [ [t6[ . 
The asymptotic form factor is predicted by perturbative 
QCD |l^,^l|. It depends on the internal pion dynamics 
that is parametrized by the quark distribution amplitude 
of the pion. The QCD expression for the pion form factor 
is |lj 

FAQ 2 )= [ dx[ dyr(y,Qy)T H (x,y,Q 2 ) 
Jo Jo 



x4>(x,Q x ), 



(12) 



where Q x = min(a;, 1 — x)Q and Q is the 4-momentum in 
Euclidean space transferred by the photon . The func- 
tion <j)(x, Q x ) is the pion wave function, that gives the 
amplitude for finding the quark or antiquark within the 
pion carrying the fractional momentum x or 1 — x, re- 
spectively. In this work we use the model for the pion 
distribution amplitude proposed by Chernyak and Zhit- 
nitsky |2^ |. This wave function was derived from QCD 
sum rules and it is written as 



</){x,Q) 



A 

2%/3 



x(l — x) 



6+ [30(2x- l) 2 - 6] 
with jj, = 500 MeV and 



a.(Q) 



72 



50 



99- 6n 



(13) 



(14) 



The other function, Th(x, y, Q 2 ), is the hard-scattering 
amplitude that is obtained by computing the quark- 
photon scattering diagram as shown in Fig. |^. The 
lowest-order expression of T^(a;,y,Q 2 ) is given by (see 
p3[ , and references therein) 



64^ (2a s [(l-x)(l-y)Q 2 } 



T H (x,y,Q) = <-- 

3Q 2 13 (1 - x){l - y) 

1 a s {xyQ 2 ) 

3 xy 



(15) 



To compute the pion form factor using the nonpertur- 
bative runningbcouplings proposed by Alkofer et al. p| 
and Cornwall Eq.(p), we solved the integrals given by 
Eq.(|l2|) substituting the quark distribution amplitude 
written in Eq.dl3|) and the expression of TH{x 7 y,Q 2 ) 
(Eq.©). 

The pion form factor result for the different forms of 
the QCD coupling in the low momentum regime is shown 
in Fig.(§). We used for Eq.(§), the lower (300 MeV) 
and the upper (700 MeV) gluon mass values for a fixed 
A = 300 MeV. These values defined the shaded area rep- 
resenting the expected range for the pion form factor, 
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^(1, y, Q) 



FIG. 2. The leading-order diagrams that contribute to the 
pion form factor. <f)(x,Q x ) is the pion wave function, that 
gives the amplitude for finding the quark or antiquark within 
the pion carrying the fractional momentum x or 1 — x. The 
photon transfers the momentum q (in Minkowski space), 
Q 2 — —q' 2 , for the qq pair of total momentum P producing a 
qq pair of final momentum P'. 

In this same figure, we also compare our results with 
the experimental data (solid line) |24j that was described 
bv the least \ 2 St (x 2 mm = 7.96742) determined in Ref. 



fit 0-46895 / _ 0.3009 \ 
Q 2 V Q 2 ) 



(16) 



The results, using the running coupling of the Eq.(||), 
agree very nicely with the experimental data for a gluon 
mass value close to 700 MeV. On the other hand, the 



calculations with Eq.(10) overestimate at least by one 
order of magnitude. 



IV. EFFECTS OF NONPERTURBATIVE 
PROPAGATORS IN THE BEHAVIOR 

In the previous section we computed F n using two dis- 
tinct forms of the nonperturbative running coupling. We 
considered that the gluon exchanged by the qq pair of 
Fig.(^) is a perturbative one. However, the SDE solu- 
tions at the same time that they give the nonperturbative 
behavior of the running coupling, they provide nonper- 
turbative expressions for the gluon propagators that at 
the origin differ drastically from the perturbative propa- 
gator. 

The large momentum behavior of these nonpertur- 
bative propagators coincide with the perturbative one, 
and, by consistency, we have to use the nonperturbative 
gluon propagators together with their respective coupling 
constants, even considering that we are computing the 
asymptotic pion form factor. So that, it is worth ask- 
ing whether our previous analysis would be distinct if we 
change the perturbative gluon propagator by the full one. 

In order to introduce this modification, we verify that 
in Eq.(|l5|) we used the perturbative QCD gluon propa- 
gator that, in the Landau gauge, is given by 



1 



Dtf) = - S . (17) 



D^)=[S„- q -^)D(q^ 



We can easily factorize D(q 2 ) in the Eq.(|T^), rewritten 
this last equation as 



Alkofer et al. 
Cornwall - m =300 MeV 




6 8 10 

Q 2 [GeV 2 ] 



Th{x,v,Q 



64tt 
~3~ 



2 -a s {K 2 )D{K 2 ) + l -a s {P 2 )D(P 2 ) 

(18) 



where K 2 = (1 - x)(l - y)Q 2 and P 2 = xyQ 2 . 

Let us now consider the two different nonperturbative 
behaviors of gluon propagators. The first one was ob- 
tained by Cornwall pM, and is given by 



Dc(q 2 ) = 



1 



q 2 + M 2 (q 2 Y 



(19) 



where M g (q 2 ) is the dynamical mass given by Eq. (CT 
The gluon propagador computed by Alkofer et al. [ |l5[ 
can be fitted by the following expression (\ 2 = 0.016) 



FIG. 3. Pion form factor computed with the different non- 
perturbative running coupling constants. The curve com- 
posed by line + square is obtained with Eq. ([l^) . The curves 
that define the shaded area are computed with Eq.(^) for the 
values of the m g — 300 (upper curve) and 700 MeV (lower 
curve). The solid line is the experimental data fit Eq,(|l6|). 
There is a nice agreement when F^ is computed with Corn- 
wall's running coupling. 



DAM 2 ) - -X^, 
<r + a 4 



(20) 



where a = 0.603 and b — 3.707. 

Once the propagators are given by Eqs. (|r^ , po| )thcn Th 
(Eq.@) will be changed to 



T H {x,y,Q ) = — - 



1 



-a s {K 2 )D A , c {K 2 ) 



-a s (P 2 )D A , c {P 2 ) 



(21) 
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We performed the integrations of Eq . ( |12[ ) numerically, 
with the amplitude Th given by Eq. (|2l| ) and their respec- 
tive running coupling constant (see Fig. ([!])). Our results 
are shown in Fig.fl). 




11 12 13 14 15 16 



CT[GeV'; 



FIG. 4. Comparison among the experimental curve (solid 
line) and our results for the pion form factor using the 
Alkofer's and Cornwall's nonperturbative propagators and 
coupling constants. 

If we compare the results of Fig.Q) with the results of 
the previous section, we can observe a striking attenua- 
tion of FiriQ 2 ) for Q 2 — > 0. It is also clear that F w (0) is 
finite for both models. This new behavior at low momen- 
tum can be understood if we notice that, for Q 2 — > 0, 



D A (Q 2 ) 
D C (Q 2 ) 





finite, 



(22) 
(23) 



in contrast with the divergent perturbative propagator. 
When we use the nonperturbative information obtained 
through the SDE with the approximations of Ref. |l5) , we 
continue to have a disagreement with the experimental 
data. In this particular case the pion form factor even go 
to zero as Q 2 — > 0. 

Obviously we should not consider this region of tran- 
ferred momentum because the kernel of Eq.(21) is valid 
only for large Q 2 , but the disagreement goes through 
the asymptotic region. On the other hand, the Corn- 
wall's propagator is still compatible with the experimen- 
tal data, but now the agreement is in favor of smaller 
gluon masses. 



V. DISCUSSION AND CONCLUSIONS 

There is an increasing phenomcnological evidence for 
the freezing of the QCD running coupling constant in the 
infrared region. It is clear that much more work has to 
be done in order to establish definitively these results. 
However, it is very satisfying to see that they are not far 



apart, and are concentrated on a region slightly below 
a s w 1. 

On the theoretical side there are many studies leading 
also to this infrared fixed point. Among these we selected 
the ones derived from the solutions of Schwinger-Dyson 
equations. 

In this work we proposed to test the compatibility be- 
tween the phenomenological values of a s (0) with the val- 
ues given by the SDE solutions. This compatibility (or 
not) can teach us if the approximations used to solve the 
SDE arc realistic or not, and if more data is accumulated 
we may even be able to discard nonphysical solutions. 

We discussed two SDE solutions for the running cou- 
pling constant and gluon propagators. One proposed in 
Ref. |Q and the other in Ref. These are the only 

ones consistent with the recent simulations on the lattice 
of the gluon propagator Jll|. These solutions have been 
obtained in Euclidean space and in pure gauge QCD, i.e. 
n f = 0. 

The effect of the number of flavors in Cornwall's so- 
lution |l4| is not so strong, and it appears in the coeffi- 
cient Pq of the coupling constant and in the gluon mass 
equation increasing the value of the frozen coupling. If a 
nonzero number of flavors produces any observable effect, 
this one should act in the same sense for both solutions. 
Therefore, we do not expect large changes in our results 
with the inclusion of fermion loops in the SDE solutions, 
and we can say that the phenomenological data on a s (0) 
is only compatible with the running coupling determined 
in Ref. @. 

It is known that different approximations in the same 
set of SDE produce different results. For example, Atkin- 
son and Bloch solved the same equations of Alkofer el. 
al using bare truncation and performing an angular aver- 
aging of the integrals. In this calculation they obtained 
a s (0) — 11.47 pq] . When the angular integrals were per- 
formed exactly they found a s (0) w 4.2 [g7|. In all these 
cases, the incompatibility with phenomenological data is 
still present. These studies can be improved requiring 
multiplicative renormalizability of gluon and ghost prop- 
agators 

It has been claimed that the asymptotic pion form fac- 
tor is quite dependent on the behavior of a s |l6[ . There- 
fore, in Section III we computed F V {Q 2 ) with both SDE 
solutions. Again, one of the solutions is clearly preferred 
than the other. Although we followed a traditional cal- 
culation performed by several authors, where the form 
factor was calculated using the nonperturbative running 
coupling, we commented in Section IV that a consistent 
treatment is obtained only if the nonperturbative gluon 
propagators are also taken into account. We modified 
the expression for the pion form factor including the full 
gluon propagator. The pion form factor is clearly modi- 
fied in the infrared in both cases compared to the result 
of the previous section. It is important to recall that the 
perturbative QCD expression for F V (Q 2 ) is not reliable 
for small Q 2 . However, for large Q 2 the incompatibility 
of one of the solutions with the data is apparent. 
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In summary, the phenomenological data on the low en- 
ergy behavior of the QCD coupling constant can be used 
to constrain the solutions of Schwinger-Dyson equations 
for the coupling constant and gluon propagators. More 
data is necessary, but the ones already existent indicate 
that some approximations made in the SDE, leading to 
a particular value of the running coupling in the infrared 
region, may not be precise enough to reveal its actual 
behavior. 
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